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Abstract. We study quantum Darwinism, the redundant recording of information 
about the preferred states of a decohering system by its environment, for an object 
illuminated by a blackbody. We calculate the quantum mutual information between the 
object and its photon environment for blackbodies that cover an arbitrary section of the 
sky. In particular, we demonstrate that more extended sources have a reduced ability 
to create redundant information about the system, in agreement with previous evidence 
that initial mixedness of an environment slows — but does not stop — the production of 
records. We also show that the qualitative results are robust for more general initial 
states of the system. 
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1. Introduction 

The theory of decoherence [HEllBj is supported by striking experimental evidence [HE [6] 
and helps explain the emergence of the classical realm in a purely quantum universe. 
Classicality is marked by several characteristics that, at first sight, appear not to be 
native to quantum mechanics. The progress over the past three decades has been due 
to the realization that these classical aspects of our Universe can arise dynamically, and 
are sensitive to the details of the systems and interactions being studied. 

Decoherence can explain why effectively classical pointer states (e.g. Gaussian 
wavepackets) are preferred for certain systems, in the sense that the system's pointer 
state is unaffected by its interaction with the environment, so that any other initial 
state will quickly become approximately diagonal in the associated pointer basis. It 
also explains why (and under what circumstances) interference between components in 
a system-environment entangled state can be ignored, the environment can be traced 
out, and the system state can be regarded as having "reduced" to an approximate 
mixture of pointer states. 

Nevertheless, there are still deep unanswered questions about the quantum-classical 
transition. In particular, decoherence alone does not explain two aspects of classical 
mechanics which we take for granted: the states of classical systems are robust and 
objective. By "robust" , we mean that observers may discover an initially unknown state 
of the classical system without disturbing it. By "objective", we mean that multiple 
observers may independently find out the state of the same classical system, that they 
will all agree on the answer, and that their measurements will leave the system in the 
preexisting (objective) state. 

In general, quantum states are neither robust nor objective even after decoherence 
has eliminated obviously quantum superpositions. When an observer measures a system 
in something other than it's pointer basis, he dramatically affects the state by re- 
preparing the system in an eigenstate of the observable he has measured. Therefore, 
when multiple observers each measure a system in different bases, they will each get 
different, incompatible results, and the last measurement will leave the system in a state 
that has little to do with the states revealed by its predecessors. What, then, explains 
the objective and robust nature of classical macroscopic objects? 

Quantum Darwinism provides the answer to this question [Til^. It recognizes that 
real observers do not typically interact directly with the system they measure. Instead, 
the system is immersed in and correlated with a (decohering) environment. The observer 
then interacts with a fraction of that environment to find out its states - to become 
correlated with the system - measuring it indirectly. For instance, when we "measure" 
the position of a chair by looking at it, our eyes do not directly interact with the chair. 
The chair's state is not affected by whether or not we open our eyes. By opening our 
eyes, we merely allow them (and hence, our neurons) to become correlated with some 
of the photons scattered by chair (and hence, its position). 

The environment in the Universe we inhabit acts as an information channel through 
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which the observer finds out about the system [9l |8]. But this information is filtered, 
as only the observables that are recorded in many copies in the environment can be 
found out from the intercepted fragment. This means that observers cannot choose any 
arbitrary basis in which to measure; they are restricted in the type of information that 
can be acquired about the system by the nature of the system-environment interaction. 
Under the condition of effective decoherence, this information can only describe the 
pointer states of the system, not superpositions thereof [I0l|9]. Indeed, the no-cloning 
theorem [TTl [T2] implies that arbitrary quantum states of the system will not be able 
to proliferate in this manner, as only certain preferred states (that turn out to be 
pointer states) can be imprinted onto many fragments of the environment [13j. In fact, 
observables complementary to the pointer states effectively become inaccessible after 
decoherence has set in since they can only be recovered through global measurements 
on the whole of the environment. 

So based only on the assumption that typical observers learn about the system 
through the environment (i.e. independent of any details about the size of the observers 
or the way they interact with the environment), we can conclude that (1) observers do 
not disturb the system ("robustness") and (2) all observers can learn only about the 
pointer basis and, consequently, will agree ("objectivity"). 

This description, which is related ^ to the formalism of quantum trajectories 
[16], can only hold when multiple observers can actually determine the state of the 
system by sampling just part of the environment. That is, information about the 
system must be recorded redundantly in the environment. This will not always be 
the true. In some cases of effective decoherence, the environment simply does not make 
multiple copies of the information (e.g. coUisional decoherence from a single high-energy 
environmental particle). In other cases, strong self-interactions scramble correlations 
and prevent information about the system from being extracted from any accessible 
fraction of the environment (e.g. coUisional decoherence from air molecules). 

The capacity of quantum Darwinism to explain the quantum-classical transition 
then rests on whether decoherence in everyday settings actually induces sufficient 
redundancy such that the classical approximations of robustness and objectivity are 
justified. This has been investigated for a spin-| particle monitored by a pure fTT] 
and mixed p!8l [T9] bath of spins and a harmonic oscillator monitored by a pure 
bath of oscillators [201 Ell- We recently showed for the first time that a physically 
realistic setting — an object illuminated by a point-source blackbody — does in fact lead 
to enormous redundancies [22]. In this work, we generalize that analysis to include 
partially- angularly-mixed illumination and arbitrary initial object wavef unctions. We 
confirm evidence in earlier studies [lU [19] that initial partial mixedness of the 
environment hampers — but does not eliminate — its ability to redundantly record the 
state of the system. We also show that quantum Darwinism in our model is robust for 
general initial states of the system. 
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2. Quantum Darwinism in Photon Collisional Decoherence 

When an object in a mesoscopic superposition is exposed to radiation, scattering photons 
wiU quickly reduce its pure, nonlocal state to a mixture of localized alternatives via 
collisional decoherence [23]. (See also [SI ESI [26l for refinements and corrections.) 
The fantastic rate of collisional decoherence has been confirmed experimentally [281 [29] . 

Observers typically access a small part of the environment (in this case, the photons 
that enter one's eye), so we will estimate how much information about the object is 
available in a subset of the environmental photons. For simplicity, we assume our 
environment consists of a large but fixed number of identical photons: S = 
where Si is the Hilbert space of a single photon in a box of volume V. We then define 
J^f = be the fragment corresponding to some fraction / of the environment 

composed of fN photons. Since each photon has the same initial conditions and 
interactions, the choice of photons with which to construct the fragment is unimportant. 
To get our final results, we will take V and to infinity while holding the physical 
photon density N/V constant. 

The primary quantity investigated will be the quantum mutual information 

^S:J=- = Hs + Hjr — Hs^j^ (1) 

between the system S and fragment J-", where H denotes the von Neumann entropy. 
From this we will calculate the redundancy Rs^ which is the number of distinct fragments 
in the environment that supply, up to an information deficit 5, the classical information 
about the state of the system. More precisely, Rs = 1/fs^ where fs is the smallest 
fragment such that ^s-.Tf = (1 ~ ^)Hs. {Hs is the maximum entropy of S. Only very 
large fragments / > 0.5 will be able to have perfect classical information, ^s-.Tf^ = Hs^ 
about the object ^IZj.) At any given time, the redundancy is the measure of objectivity; 
it counts the number of observers who could each independently determine the state of 
the system (up to a small residual uncertainty 6) by interacting with disjoint fragments 
of the environment. 

Following Joos and Zeh [23j , we take our system to be a dielectric sphere of radius 
a and relative permittivity e in a pure state. The system and the photons in the 
environment are assumed to be initially unentangled: p° = Pg • • • Pg, where 
ps and pe are the density matrices of the system and of a single photon, respectively, 
and a superscript "0" denotes prescattering states. 

The system is illuminated by photons originating from a far away blackbody of 
temperature T that covers B C S, where S is the unit sphere "sky" as seen from S. Let 
Q < An he the solid angle measure of B. See figure [ij 

In the Hilbert space of a single photon, we break the momentum eigenstates into 
a tensor product \k) = \ki)\h)/k of magnitude and directional eigenstates. (The factor 
of 1/k comes from the Jacobian determinant associated with fact that the kets in this 
infinite-dimensional Hilbert space are densities, not normalized vectors.) Blackbody 
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Figure 1. A dielecric sphere of radius a and permittivity e is initially in a superposition 
with separation Ax = \xi—X2\. The object is subjected to radiation from a blackbody 
at temperature T that originates from a patch B C S (where S is the unit sphere "sky" 
as seen from S) with solid angle Q < Att. The complement of B is B. The photons 
propagate in directions labeled by n, which makes an angle 6 with the vector Ax. 



illumination is then described by 

pI = / dkp{k)\k){k\® / — (2) 

where p{k) oc k'^ /[exp{kc/kBT) — 1] and c is the speed of light. Above, the probability 
distribution over the solid angle is the normalized characteristic function of B. This 
"step- function" illumination is not contrived. Perfect blackbodies are Lambertian, which 
means that surface elements appear to have the same brightness no matter the angle of 
viewing. In other words, the sun appears as a uniform disk of illumination in the sky; 
it is not dimmer near the edge. 

We ignore the self-Hamiltonian of the object and assume it is heavy enough to have 
negligible recoil from photon scattering, so the evolution is governed by 

Ut\x)\k) = \x)Ss\k), (3) 

where is a scattering matrix acting on the single photon state when the particle is 
located at x. Elastic scattering leads to 

Ss\k) = Ss\k)\h)/k = \k)Sl\n)/k, (4) 

under the magnitude-direction decomposition of the photon momentum states. 



2.1. Decoherence 

For now, we take the initial state of the object to be a (Schrodinger) "cat" state: 
Ps = \^){^\^ \^) = (l^i) + I ^2)) /a/2- The decoherence of the superposition is governed 
by the decay of the off-diagonal terms in the position basis, 

\{xi\ps\x2)f = 7^ \{xi\ps\x2)\^ , (5) 

where 

2 

(6) 



7 



Tr 
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The complex value 7 is the decoherence factor attributable to a single photon. The 
two-dimensional ps can be diagonalized and its entropy is 

//5=ln2-V— (7) 

n=l ^ ^ 



= In 2 - a/T arctanh Vf - In Vl - T (8) 

where F = 7^ is the decoherence factor associated with the environment as a whole. 
The function 



00 



h(x) = \fx arctanh \fx + In \J\ — x = — (9) 

2n(2n — 1) 

will appear often. Note that /i(0) = 0, h{l) = In 2, and 

x/2 < h{x) < xln2. (10) 

Also, h{x) is analytic and monotonic on the interval [0, 1] and, for small x, h{x) ^ x/2. 
Applying ^ gives 



Tr 



S^^yeSl]= dkp{k) J^^ (11) 

To get the key matrix element appearing on the r.h.s., we use the classical cross section 
of a dielectric sphere [30] in the dipole approximation (A ^ a) and assume the photons 
are not sufficiently energetic to resolve the superposition individually (A ^ Ax) [fj This 
gives [22] 

(-1^241-) = 1 - ^^(3 + llcos^^)^^^^ + O (y-) , (12) 

where a = a[{e — l)/(e — 2)]-^/^ is the effective radius of the object, 9 is the angle between 
h and Ax, and t is the elapsed time. 

For increasing photon momentum eigenstates become diffuse so individual 
photons decohere the state less and less. In other words, 7^1 because I. This 

is balanced by an increasing number of photons in the box. In the limit ]/, ^ oo, we 
combine ([g]), (11), and (12) and use e = limg^oo(l + 1/?)^ to get the decoherence factor 

F= lim 7^ = e-'/^^, (13) 
where td is the decoherence time. It's inverse is the decoherence raie[§} 

--'4(l)(li)^/.f(— (") 



(15) 



I High energy photons are uninteresting as they can easily "see" the separation Ax and each will 
become completely correlated with the position of the object; the redundancy will be just be equal to 
the number of scattered photons. 

§ This is related to Joos and Zeh's convention by r^^ = 2r~^ = 2AAx^, where r and A are the 
"characteristic time" and "localization rate" of [23]. 
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where ({n) denotes the Riemann zeta function, which arises from integrating over the 
thermal distribution. In the second hne, we have expressed the number density of 
blackbody radiation in terms of the apparent sohd angle, N/V = J7("(3)/c^T/(27r^c^/i^), 
and we have rewritten the integral to emphasize that it is merely the average of the 
trigonometric quantity 3 + llcos^^ over B. If we ignore the order-unity change in 
that average (it is obviously constrained to lie between 3 and 14) then the decoherence 
rate is linear with the apparent size of the blackbody. This is natural because thermal 
radiation is uncorrelated. Each new photon contributes an independent multiplicative 
decoherence factor, which combine additively in the decoherence rate. 

For blackbodies that are far enough away to be approximated as point sources, the 
irradiance I (radiative power per unit area) is a more physically accessible quantity 
than the the solid angle f], especially in the presence of optical distortion. In that case 
we use N/V = (//c/c^T)[2!C(3)]/[3!C(4)], to get the point-source decoherence rate [22] 



. /47r\ /8!C(9)\ , ,^Ja^Ax^k%T^ , , 

where 9 is the angle between Ax and the point source. 

At the opposite extreme where Q = 47r,B = §, we recover the decoherence rate for 
isotropic thermal illumination [26] : 

_ . / 8!C(9)\ a^Ax^k%T^ 

This is the decoherence rate when the system is surrounded by a uniform blackbody, 
e.g. inside an oven. For the corresponding decoherence time, we have introduced the 
symbol = ti)|b=§ since it will serve as a useful B-independent timescale in the rest 
of this paper. 

For concreteness, consider a blackbody that appears as a disk on the sky centered 
on a direction z which makes an angle x with Ax. That is, B = {(0, 0) G S | ^ < ^o} 
for some maximum polar angle ^o- See figure [4] For such disks, the decoherence rate is 

T^^ = — [40- cos^o(51 - 33cos2x) + cos^^o(ll - 33cos^x)] T-\ (18) 
80 

This is plotted in figure [4](a) in terms of the solid angle Q = 27r(l — cos6^o)- It is 
monotonically increasing with Q since each additional photon can only further decohere 
the system. 



2.2. Quantum Darwinism 

To get the redundancy in the environment, we will need to find the mutual information 
Is:T = Hs + Hjr — HsT- We can avoid calculating Hst by using the identity [Eq. (8) 
of [19]] 

Is:T = [Hj- - H^] + [HsdE - HsdE/j] , (19) 
where Hsds = Hs is the entropy of the system as decohered by the entire environment 
^, and Hsds/T is the entropy of the system if it were decohered by only S/J^. We get 
Hsds/T from Hs^ ([7|), by making the replacement F r^~f . 
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The calculation of Hjr is tedious, and we relegate the details to the appendix. The 
change in the entropy of the fragment is 

^Hf = Hf-H% = \xi2-Y^ -, (20) 

^ ^ :f Z^2m2m-1)' ^ ' 

where the effect of the initial mixedness of the environment on the production of records 
is accounted for by the single parameter 



/b dn4dm|^( 



n, m 



12 



/b dn/g dm|^(n,m)|2 

For reasons that we explain below, we call a the receptivity of the environment with 
respect to the decoherence process. Above, B = S\B is the complement of B and 
we define \g{n^m)f = J dkp{k) \{n\{S^^^ S^^ — I)\7fi)\ . The function \g{h^m)f is a 
measure of the distinguishability of the out states S^^\h) and S^^\m) for different 
incoming angles h and m of £ and for different locations Xi and X2 of 5. The receptivity 
a is a dimensionless ratio constructed from this function that, from the form of (21), 
we know obeys < a < 1. 

Now that we have the change in fragment entropy AH^^ = Hp — i/^, we use (19) 
to finally write down the mutual information 

^ p(l— /)m pa/m pm 

Xs-^ = In 2 + V ^ . (22) 

^ 2m(2m-l) ^ ^ 

m=l ^ ^ 



The summations in (22) can be written in a closed form analogous to ([8|), but the power 
series is more useful for calculating the redundancy. For large times, F = exp(— t/r^) is 
exponentially small and the sum is dominated by the lowest power of F. If < / < 1/2 
and a 7^ 0, then < a/ < (1 — /) < 1 and 

Is:^^ ^ ln2-^r"^. (23) 

So long as the information deficit is not unreasonably large, 6 < 1/(2 In 2) ^ 0.72, we 
can estimate the redundancy in the limit t ^ td: 

where 

T^' = aT^' (25) 

is the redundancy rate — the characteristic rate at which records about the state of the 
system are produced 

The redundancy ([24]) depends only weakly (logarithmically) on the information 
deficit 5^ which is consistent with previous results [191 EI]. The redundancy increases 



II The times for which this is a good approximation to the true redundancy are shown in figure [s] For 
added rigor, we can use (10) to get Xs-.jr > ln2(l — T^^ — F), which, for t > rdln{2/6), yields this 
lower bound on the redundancy: Rs > (t/r/))/ ln[((5 — F)"^]. Since F decays exponentially in time, this 
conservative bound tracks our estimate (24) very closely for all S < 1/(2 In 2). 



Redundant Information from Thermal Illumination 



9 




Figure 2. The quantum mutual information ^s-.j^f (22) versus fragment size / at 
different elapsed times for an object illuminated by blackbody radiation. Individual 
curves are labeled by the time t in units of the decoherence time td (14). The same 
time slicing is used for all values of the receptivity a. Top left: a = 1. For t <td (red 
dashed lines), the information about the system available in the environment is low. 
The linearity in / means each piece of the environment contains new, independent 
information. For t > td (blue solid lines), the plateau shape of the curve indicates 
redundancy; the first few pieces of the environment reveal a lot of information about 
the system, but additional pieces just confirm what is already known. On the plateau, 
the mutual information approaches it's maximum classical value, Hs = 1 bit = In 2 
nats. The remaining information (i.e., above the plateau) is highly encoded in the 
global state, in the sense that it can only be read by capturing almost all of £. Top 
right: a = 0.5. Redundant copies are only produced at half the rate and the mutual 
information is no longer anti-symmetric about / = 0.5 because the environment is 
mixed in the information-storing degrees of freedom. Nevertheless, the plot approaches 
the same symmetric plateau shape for large t, illustrating that extensive redundancy 
is still achieved. Bottom left: a = 0.01. For receptivity this low (which is only 
expected for nearly isotropic illumination) the mutual information is initially greatly 
skewed. Still, this just slows down acquisition of information by a factor of = 100. 
(Compare the t = lO^-^ri^ curve for a = 0.01 to the t = lO^-^ri^ for a = 1.) 
Bottom right: a = 0. Only for this idealized case of perfectly uniform illumination is 
information storage halted. This is because the directional photon states are already 
"full" and cannot store more information about the state of the object. Zero redundant 
copies are produced and the mutual information approaches as t ^ oo for all / < 1. 



linearly with time at a rate proportional to the decoherence rate. This is intuitive 
because (1) photons scatter off the object at a constant rate and (2) it is precisely the 
dependence of photon out states on the position of the object (roughly corresponding 
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Figure 3. The redundancy Rs versus time for an object illuminated by blackbody 
radiation of varying apparent solid angle, with information deficit S = 0.01. The exact 
redundancy (solid black lines) is plotted starting when Rs = i.e. when it first become 
possible for an observer to determine the state of the system to within the information 
deficit 6 by accessing the entire environment. For times t ^ td^ the redundancy is 



well approximated by the linear estimate (dashed blue line) in (24). 



to a record) that causes decoherence. 

It can now be seen why we caU a the receptivity of the environment; it determines, 
for a fixed rate of decoherence, the rate at which records about the systems are created 
in the environment. For maximum receptivity (a = 1), the redundancy rate is equal to 
the decoherence rate and records are produced at the maximum speed. For vanishing 
receptivity (a = 0), no redundant copies are produced no matter how effective the 
decoherence. 



2.3. Receptivity 

We now examine how the receptivity depends on the distribution of the iUumination and, 
through the operator S^_^^ S^^^ the scattering matrix. However, most of our discussion 
wiU not rely on the particular angular dependence of the differential cross section. 

In general, both the decoherence rate r^^ and the receptivity a will vary with 
different choices of B. The decoherence rate r^^ changes for two reasons: (a) Different 
incoming photons contribute independently (that is, additively) to the decoherence rate, 
so larger or closer blackbodies covering more of the sky will naturally decohere the 
system faster, (b) Less dramatically, the contribution per unit solid angle to the rate 
of decoherence changes as the blackbody is moved around the unit sphere because of 



the factor of 3 + llcos^^ in (12). The receptivity a similarly depends both on (a) the 



solid angle of the regions integrated over in (21), and (b) the angular dependence of the 
integrand \g{h^rh)\^. 
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To disentangle these two quantities, simply consider a blackbody dimmed by some 
uniform intermediate medium. Let /3 G [0, 1] be the fractional degree to which the 
intensity of the illumination is reduced, expressed by modifying the number density 
N/V /3N/V (equivalently, the intensity / /3I). The reduces the decoherence 
rate accordingly, r^^ I^^d^ ) but leaves the receptivity a unchanged. We can then 
consider an arbitrary family of blackbodies B^^^ with different a^'^\ taking i = to be 
the one with the largest decoherence rate: 4°)-^ > T^S~^ for all i. If we dim each 
B^^^ by = ^ we equalize all the decoherence rates without changing 

the receptivities. This gives a physical interpretation for considering, with a fixed 
decoherence rate, how the receptivity (and hence the redundancy rate) depends on 
the shape of the blackbody. 

In the appendix we show that \g{n^ which is bounded, has angular dependence 

\g{n,m)\^ oc (1 + cos^ (9^,^)(cos (9Ax,n - cos (9Ax,m)^ (26) 

where cos 0^,6 = a • 6 is the cosine of the angle between the unit vectors a and h. The 
explicit form of the receptivity for general B is then 

_ /jg dn dm (1 + COS^ (9^^^) (cos 6>Ax,n - cos ^Ax^m)^ . . 

Jjg dn Jg dm (1 + COS^ On^rn) (cos 6>Ax,n - COS 6>Ax,m) 

We know that < a < 1 and, since \g{n^ m)p only vanishes when cos ^Ax,n = — cos ^Aa:,m 
(a set of measure zero on § x §), we see that the two extremes are realized only in the 
following physical situations: 

• For B = {no} (point source illumination from the direction fio), = 0, a = 1, and 
we recover (13) of 



Ep(l-/)m p/m 
-r ^ • (28) 
2m 2m - 1 ^ ^ 

See figure [2} Of course as noted above, exact point sources are a mathematical 
idealization; total fiux of blackbody radiation is proportional to fi, so this is really 
the case when finite-size eflFects can be ignored. 

For B = S (isotropic illumination), = 47r, a = 0, and we recover (17) of [22]: 

'^s-.j^f = Hsds — Hsds/T (29) 

= L 2n(2n - 1) ' 

n=l ^ ^ 

For times t ^ r^^ the entropy of the system as decohered by the entire environment, 
Hsds^ and by just the complement of the fragment, Hsds/T^ are both exponentially 
close to In 2, so that Ts-.Tf = Hsds — Hsds/T is only non-negligible for the brief period 
when t ^ td- This is plotted in figure [2| which shows that the mutual information 
barely rises from zero before fading away, never yielding a single redundant copy. 



(For strictly vanishing a, the approximation used for (24) breaks down.) The 



photon directional states, which are the component of the environment in which 
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information about the object is stored, are initially fully mixed and so cannot hold 
any new information. In other words, an observer relying on scattered radiation 
in an oven can see nothing [32j. This makes it clear that decoherence — which is 
maximized for isotropic illumination — is not sufficient to guarantee redundancy [8] . 



Between these two extremes, the form ( |21| ) indicates that the receptivity a will 
tend to decrease with increasingly large B In physical situations (e.g., objects lit 
by light bulbs, the Sun, or ambient light), we expect illumination to be nonuniform. 
This will correspond to receptivity that is not particularly close to either or 1, so 
that the initial mixedness of the photon environment decreases the redundancy only 
by roughly a factor of order unity (in accordance with detailed calculations made of 
spin- 1 systems [E]). Since even very tiny objects have extremely short decoherence 
times [Il[23], redundancies will still be very large for realistic illumination. 

To get intuition about the dependency of a on B, we again specialize to the case of 



the blackbody disk with solid angle Q so that the integrals in (27) can be preformed. (See 
the appendix.) The results for the decoherence rate, the receptivity, and redundancy 
rate have been plotted in figure [4] 

The receptivity of the decoherence is related to the haziness ^iHl [H] (initial entropy) 
of the environment. Environments with zero haziness in the information-storing degrees 
of freedom (photon directional states) will have maximum receptivity. Environments 
with maximum haziness will have zero receptivity. However, the receptivity is not a strict 
function of the haziness since it depends also on the form of the scattering operator, i.e. 
the way in which S makes its mark on S. 

2.4' Origin of difference between decoherence and redundancy rates 

The key difference between the redundancy rate r^^ and the decoherence rate r^^ is the 
type of matrix elements on which each depend. Mere decoherence of the two pointer 
states \xi) and \x2) of S is determined only by the overlap of the 5-conditioned states 
of 5, p^^^^ and p^^^^ Initial mixedness of the environment is accounted for by simply 
averaging the overlap of the 5-conditioned pure states of f , S^^\n)^ when computing 
the decoherence factor: 

2 



7 



Tr 



dkp{k)J^f{n\Sl^S:^ 



n) 



(31) 



So long as this average overlap is unchanged, mixedness of the environment does not 
affect decoherence. Importantly, the decoherence rate depends only on inner products 

^ Even for a strictly enlarging sequence of blackbodies M{s) (where B(5i) C B(52) and Vt{si) < ^{32) for 
si < S2) the receptivity does not necessarily decrease monotonically with For the dipole scattering 
cross section considered in this work, as well as a few other example scattering operators, we have been 
able to construct unusual examples with Bi C B2 and < Q2 such that ai < a2. However, simple 
blackbody shapes like uniform disks have receptivity that decreases monotonically with solid angle Q. 
The number of local extrema a takes with growing B is restricted by the size of higher frequency terms 
in the Fourier expansion of |^(n, m)p and the complexity of the shape of B. Trivially, if |^(n, m)p were 
constant, then the receptivity would just decrease linearly with Q. 



Redundant Information from Thermal Illumination 



13 




Figure 4. Decoherence and quantum Darwinism for an object illuminated by a 
blackbody disk. Left: The disk has solid angle Q and the center of the disk z makes 
an angle x with Ax. Precisely, B = {{0^(j)) e E> \ < Oq} where cos^o = 1 — ^/27r. 
Right: The decoherence rate, the receptivity, and the redundancy rate for a blackbody 
disk as a function of disk size. We consider the three cases X = 0° (green solid line), 
45° (blue short-dashed line), and 90° (red long-dashed line), (a) The decoherence rate 
r^^ (18) normalized by the maximum decoherence rate T^^ obtained ai ft = An. The 
decoherence rate is monotonically increasing with a property that holds regardless 
of the scattering cross section, (b) The receptivity a (A.41) decreases monotonically 
with Q (i.e. with increased angular mixing of the environment), a relationship that 
holds for arbitrary x- (c) The redundancy rate = <^r^^ normalized by T^^. The 
redundancy rate is symmetric about Q = 27t because B and its complement B appear 
on equal footing in ( A.31 ). For = 0, there is zero illumination and so there is neither 
decoherence nor quantum Darwinism. For Q = 47r, the decoherence is maximized 
but there is no Darwinism because the receptivity vanishes; the angular mixing of the 
environment is total, allowing no recording of information. Away from either extreme 
the redundancy rate is within an order of magnitude of the decoherence rate. 
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\m 



X2 



Figure 5. The various out states for different initial states of £ and different scattering 
states of S. Let n,m G B be possible directions for incoming photons and let Xi and 
X2 be the pointer states of the object. (We ignore photon momentum, which decouples 
since the scattering is elastic.) Then there are four distinct out states and therefore 
four non-trivial inner products. The decoherence rate r^^ [see (A.30)] depends only 
on inner products between like initial states of the environment, {h\S^^^ S^Jh) and 
{7fi\S^_^^ S^Jifi), but the redundancy rate [see (A. 31)] depends also on {h\S^_^^ S^Jm) 
and {h\S^J S^^Itti). High redundancy depends not just on small overlap between the 
<S-conditioned states for the range of initial states of £ but also on the observer's ability 
to infer the state of S from imprints left on different initial states of £. 



between the 5ame initial states of S conditioned on different states of S. 

On the other hand, the production of records is very sensitive to the initial 



mixedness. The redundancy rate is proportional to the numerator of (A.26). That 
numerator includes matrix elements of the form | (nlS'l^^S'l^ |m) | where, crucially, 
h 771. For redundant records to be produced, there must be small overlap between 
5-conditioned states for different pure states of £ in the initial mixture. In other words, 
the observer must be able to distinguish the imprint of the pointer states of the system 
on different initial environment states. When [(nl/S^^'^'S'l^lm)! is not small for n 7^ m, 
then the observer cannot tell whether she has sampled (a) a photon that started in state 
|n) and scattered off the system in state \xi) or (b) a photon that started in state \m) 
and scattered off the system in state \x2). See figure [5| 

3. General superpositions 

So far we have considered only objects localized in a balanced "cat" state: ^ 
[6{x — Xi) + 6{x — X2)]/2. We now gather some evidence that relaxing this assumption 
to allow more general initial object states does not give qualitatively new behavior. 

For an arbitrary initial wavefunction the decoherence behavior is very simple. 

There is no self-evolution of the object, and the initial off-diagonal terms ^^(x, x') decay 
exponentially in time at a rate that is a function of Ax = |x — x^|. As described in [25j, 
the rate is proportional to Ax^ for small distances and saturates to a constant value 
for large distances. We would like to know if the mutual information Is-.j- is similarly 



+ The apparent dependence of (A.30) on matrix elements of the form {n\S^_^^ S^Jrh) for n 7^ m can be 



removed by using the completeness relation Jg dm|m)(m| = /. The same cannot be done for (A.31). 
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well-behaved for general initial object wavefunctions. 

First, we consider the mutual information for the "unbalanced cat" state 

illuminated by a point-source. The new post-scattering density matrices are 



a,b=l,2 



N 



and 



fN 



(32) 

(33) 
(34) 

(35) 



a=l,2 



Diagonalizing is not any more difficult than the balanced case, yielding the eigenvalues 
1 



X2 



2N 



= i±^V/^ + (i-/i)r, 



(36) 
(37) 



where /x = (pi — ^2)^. From this, we see i/^ for the unbalanced cat can be obtained 
from the balanced case by making the replacement F ^ /x+ (1 — /x)F. (The decoherence 
factor, defined as the inner product between the environment states conditioned on pure 
system pointer states, is still F.) 

We decompose the photon Hilbert space as in the appendix, (A.l), only now with 

fN 



P 



(38) 



1=1 



i=l 



It is only a little more work to show that the eigenvalues of (which carry momentum 
dependence) are similarly modified as 



(39) 



Carrying out the momentum integrals leads to the similar replacement F-^ ^ /i + (1 — 

The net effect of unbalancing the cat state on the mutual information is to make the 
replacement x /x+(l— for x = F, F-^, or F-^"-^, in (28). The only qualitative change 
to the partial information plot is to lower the classical plateau to the new maximum 
system entropy, Hs = —pilnpi — p2lnp2 < In 2, and to "soften" the shoulders. (See 
figure |6|) The mutual information for the extreme case (/i = 0) is trivially zero, but the 
limit of the renormalized mutual information for extremal probabilities is 



lim 

12^1 



HZ 



= 1 + _r/ -r. 



(40) 



Note that there is finite softening as /x ^ 1. 
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Figure 6. The partial information plot of unbalanced cat states, for t = 1{)td and 
II = (pi — P2)^ = (long-dashed blue line), /i = 0.5 (medium- dashed green line), 
/i = 0.9 (short-dashed red line), and /i ^ 1 (solid black line). More unbalancing 
(increased fi) moves the classical plateau down to the new maximum system entropy, 
Hs = —pilnpi — p2lnp2, and softens the shoulders — but with only finite softening 
as /i ^ 1. Left panel: unnormalized mutual information. (The /i ^ 1 case is not 
visible.) Right panel: renormalized mutual information, to clearly show softening of 
the shoulders. The inset shows the renormalized t = 100 td case. 
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Figure 7. The partial information plot of balanced M-way cat states for t = IOtd and 
M = 2 (long-dashed blue line), M = 3 (medium-dashed green line), M = 10 (short- 
dashed red line), and M = oo (solid black line). Increasing M moves the classical 
plateau to Hs = InM and softens the shoulders, but with only finite softening as 
M ^ oo. Left panel: unnormalized mutual information. (The M = oo case is not 
visible.) Right panel: renormalized mutual information, to clearly show softening of 
the shoulders. The inset shows the renormalized t = 100 rjj case. 



Now we look at a "balanced M-way-cat" state, ip^x) ^""^^^ X^fli ^(^ ~ ^a), 
exposed to a point-source blackbody. Let us assume the special case where the 
decoherence factors for all off-diagonal elements of ps are equal. That is, assume 
{h\S^J S^^\h) is the same for all a 7^ 6 (such as when M = 3 and Xi, X2, and Xs 
form an equilateral triangle in a plane perpendicular to the direction of illumination). 
Generalizing from the original case of the balanced 2-way-cat does not require any new 
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tricks. The corresponding mutual information is 



which can be checked to reduce to ( |28| ) for M = 2. 

This yields a partial information plot with slightly softer shoulders and, as expected, 
a classical plateau at Hs = InM. (See figure [7|) The limit of the renormalized mutual 
information for large M is 

lim ( 5^^) = 1 + Vf'~^ _ ^ _ . (42) 



Hs 



Again, there is only finite softening for even the extreme case (M oc). 

Balanced M- way-cat states with arbitrary decoherence factors are difficult to handle 
because they require diagonalizing a matrix with M[M — l)/2 different off-diagonal 
terms. We can still say the following. Let 7^^, for i ^ j, be the set of decoherence 
factors. Call the factor with the largest (smallest) absolute value 7w (ts), with "W" 
("S") standing for weak (strong) decoherence. Let (p|) be the hypothetical state 
resulting from setting all decoherence factors to 7w (ts)^ aiid likewise for p^jr {p%jr) and 

There is strong numerical evidence that we can bound H'g < Hs < Hg, in 
agreement with intuition. The same is true for Hjr and Hst because, in the proper 
basis, pjf and p^jr take the same form as ps- Finally, we note that in the large-t (small- 
r) limit, 

0<Hl-Hj =0(rw), (43) 
0<H^j.-Hj =0(r^), (44) 
0< Hl^-Hf-, = 0{T'^f). (45) 

This limit is quickly reached since F = e"^/'^^, and this is exactly where we would like 
to calculate redundancy. For times t long enough to make higher order powers of F 
negligible, the mutual information Is-.Tf — and hence the redundancy Rs — is bounded 
by the associated values for the strong and weak decoherence: 

lIr,<^S:T,<llr,, (46) 
Rf <Rs< Rl (47) 

for /< 0.5. In other words, having unequal decoherence factors for a balanced M-way 
cat state does not drastically alter the behavior of the mutual information; for times 
much larger than the decoherence time r^, the redundancy can be bounded from above 
and below by calculations using the largest and smallest decoherence factor. 



For "unbalanced M-way-cat" states, ip{x) = '}2^=iPa^{^ ~ ^a), (43 - 45) still 



hold, but we cannot use (41) to get a power series for the entropies in terms of F. 
But we do know that the classical plateau must still form at Hs = —'}2^=iPa^^Pa 
since, for large times, all off-diagonal terms are driven to zero. As M ^ cx), these 
states can approximate a generic continuous wavefunction, but subtleties enter. The 
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maximum entropy Hs = InM diverges, while the decoherence factor associated with 
two sufficiently nearby points becomes arbitrarily weak. For these small distances, 
the self-Hamiltonian of the system will no longer be negligible on the times scales of 
decoherence. 



4. Conclusion 



We illustrate quantum Darwinism [3 [9l El |33] in systems decohered by blackbody 
radiation. Such radiation is by far the dominant form of illumination in everyday life and 
is the medium though which we, as observers, gather most of our information. The huge 
redundancy growth rates we have calculated support the claim that a purely quantum 
universe can account for the appearance of objective and robust classical states. 

This is the first realistic model of quantum Darwinism, and the first with an 
environment with the capacity for two distinct types of mixing. The most important type 
of mixing is in that component of the environment responsible for storing information 
about the system (the angular degrees of freedom). We have shown, in agreement with 
previous studies of abstract systems [iHl [IS], that mixing of this type decreases the 
environment's ability to record information about the system — without decreasing its 
ability to decohere. The other type of mixing (the energy spectrum) only effects records 
insofar as certain modes (higher energy) are better able to resolve, and therefore record, 
states of the system. 

We have extended the model to more general initial states, giving evidence that the 
qualitative features of the mutual information and redundancy are robust. 
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Appendix 

In this appendix we calculate Hjr (20), the entropy of a fragment of the photon 



environment S. We take advantage of the special form of our model; the elastic scattering 
conserves energy but mixes photon direction. This alows pjr to be decomposed into k- 
blocks: 

P^ = J dXP(x) ® pj , (A.l) 

where x = (^i, • • • , ^/a^) is the vector of the magnitudes of the photon momenta 
of J^, p{x) = Yli=iP{^i) the momentum spectrum probability distribution, and 
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= The normalized p^, which hves in the Hilbert space of 

angular eigenstates, is the density matrix conditional on the particular set of momenta 
X- 

Since pj=- is block-diagonal, it's entropy is 

= fNH\p{k)] + J dxpix) , (A.2) 

where is the entropy of and II[p{k)] = II[p{x)]/ f N is the entropy associated 
with the energy distribution p{k) of a single thermal photon (which diverges since the 
photon Hilbert space is infinite dimensional). 

Handling the mutual information is easier if we discretize the angular directions. 
Let Af] be the solid angle associated with a single discretized state and = Q/AQ 
be the dimension of the projector onto the directions in B. Then send 



J VAn 



\n). (A.3) 



In the discrete picture, the conditional angular states are 



X 1 
pi = - 



P±Q 



(A.4) 



where 



Q = ®Q'^\ (A.5) 
P« = J]5j;|n),(n|5j;t, (A.6) 

Q^^ = Y.^%\n)MS%' . (A.7) 

nGB 

Notice that P^'^^ and Q*^*-' are DB-dimensional projectors unitarily related by [7*^*^ = 
S^^^S^^y P and Q are likewise unitarily equivalent and therefore they can be written in 
an appropriate basis as [34j 

where C and S are commuting positive matrices obeying C'^+S'^ = I. (Their eigenvalues 
are the cosines and sines of the canonical angles between the subspaces into which P 
and Q project.) The eigenvalues of = (P + Q)/2D^ are given by (1 ± \C\) /2D^ . 
To calculate the C, it turns out that it is easier to diagonalize the related matrix 

PQP = (g) p«g«p« = (g) ) = ( "o' ) • ^^-'^ 

* There is a technical comphcation that subspaces on which P and Q commute (which necessarily exist 
when Dm is odd) must be handled separately. This is straightforward, with no changes to the final 
results. 



Redundant Information from Thermal Illumination 



20 



By (A. 6) and (A. 8), the first subspace corresponds to the span of the basis {|n) = 
Si^\n) I n G B}. The elements of P^'^) Q^'^) P^'^) in that subspace are given by 



= ^^^{n\m){m\n') + {n\B\n') 



(A.IO) 

(A.ll) 
(A.12) 



By (12), the elements of the Hermitian matrix B defined above are of order V ^. Call 
the eigenvalues h{j)^ with j = 1, . . . , D^. The eigenvalues of are then 



2D 



fN 



1 ± n + 



(A.13) 



which are indexed by the vector J = {ji, . . . ,jfN) and the sign ±. When we plug these 
into the entropy formula and Taylor expand the logarithms, we get 



In 2 



9 DfN XI XI 



In 



1 ± JJ + Kji 



In 2 



D 



ln2- 



D 



2m(2m - 1) 
2m(2m- 1) Ai 



E 



m 
s 



(A.14) 
(A.15) 

(A.16) 
(A.17) 



Because the b(j) are of order V, we need only keep the s = 0, 1 terms in the sum over 
s. It then makes sense to define 



1 
1 



Tr 



p{i)Q{i)p{i) _ |^)(^| 



met 



We now integrate over x to get the full entropy for p-p: 



AH 



oo ^ . 



2m{2m - 1) 



oo 



i + ^E^(^'^)+^(^" 



/»oo 

l + m / d/i;p(A;)^(fc) + 0(y-2) 

^0 



(A.18) 

(A.19) 
(A.20) 

(A.21) 
(A.22) 

(A.23) 
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oo ^ 
m=l ^ ^ 

with Z = dkp{k)z{k). Then we define 

a = (A.25) 
in 7 



(A.26) 



so that 



exp(mZ/A^) = r*""-^ = 6-"^"^*/^°. (A.27) 
Moving to the N,V ^ oo hmit, we get 

A/i"^. = - M = In 2 - V ^ (A.28) 

^ ^2m2m-l) ^ ^ 



and 



12 



/b dn/g dm|^(n,m)P' 

where \g{n^m)f = J dkp{k) \ {n\{S^^^ S^^ — /)|m)|^. Note that |5'(n,m)p is symmetric 
in h and m for any scattering operator that is rotationally invariant. 
In passing, we note that, from ([2T|), 

and 



^ (X dh dm|^(n,m)|^ (A.30) 
Jm Je> 



.-oc/"«/d»M„,™)|^ (A.31) 

Allowing bars to denote quantities associated with the complementary blackbody B, 
this means that the two decoherence rates sum to that of isotropic illumination, 

T-' + f^'=T-\ (A.32) 

and that the redundancy rate is invariant when interchanging B ^ B , 

= oi/td = OijfB = (A. 33) 

Previous studies of collisional decoherence [23l EU [25l [26l [27j required the 
calculation of l^h\S^_^ S^\(m) only for h = tti. We will use most of the same tricks 
to handle the case of n 7^ m: 

14^-51.1^)1' = \{n\{I + .)e'^--^^/\l + + • • ■)\m)\^ (A.34) 
= \{h\Tk\m)\^ + \{h\Tk\m)\^ (A.35) 
- (^{h\Tk\m){m\Tk\h)e'^''-^^"^'^ + c.c.) + • • • 

= 2\{h\Tk\rn)\^ 1 - cos (^Ax ■ k{n - m) / +•••, (A.36) 



(n 
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so 



\g{h,m)\'= I dkp{k)\{h\Sl^^Sljm)\' (A.37) 



dkp{k) 



(l + cos^ On^m) ' k{n — rh)/h) 



+ 



(A.38) 



OC (l + COS^ On^rn) (cOS (9Ax,n " COS <9Ax,m)^ H , (A. 39) 

Above, T/e is the Hermitian operator that generates the unitary S^^^ for scattering off an 
object located at the origin, x = 0, and P is the momentum direction operator. EUipses 
denote higher-order terms ignored in the dipole (A ^ a) and short-separation (A ^ Ax) 
hmits. In the last line we have isolated the angular dependence since all other factors 
will cancel in the formula, ( |A.29 ), for a\ 

^ _ dn dm (1 + cos^ (9^^^) (cos Oax^ti - cos 9Ax,mf 

dh /g dm (1 + C0S2 9n,m) (cos 9Ax,n - COS 9Ax,mf 

This formula allows one to calculate the exact receptivity for any blackbody B for which 
the integrals can be performed. 

If we specialize to the case of a blackbody disk, B = {(^, 0) G § | ^ < ^o}, whose 
center makes an angle x with Ax, then 

a = [-117c^ + 2954 - 575c^ + 685 + GcJ (21c^ - 55^ + 135c^ + 75)] x (A.41) 
{ce + 1) [32 (40 + llc,(l + c,) (3^ - l))] 

where we have abbreviated cq = cos^o, = cos^. This can be combined with 
the previously calculated decoherence rate (18) to get redundancy rate r^^ = c^r^^ 
for the disk. All three quantities are plotted in figure [4] in terms of the solid angle 
f^ = 27r(l-cos^o). 
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